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Abstract 

o ■ 

CN| ' We use interacting particle systems to investigate survival and extinction of a 

$_i . species with colonies located on each site of Z d . In each of the four models studied, an 

individual in a local population can reproduce, die or migrate to neighbouring sites. 
. We prove that an increase of the death rate when the local population density is small 

OO ' (the Allee effect) may be critical for survival, and that the migration of large flocks of 

t— I . individuals is a possible solution to avoid extinction when the Allee effect is strong. 

We use attractiveness and comparison with oriented percolation, either to prove the 
extinction of the species, or to construct non-trivial invariant measures for each model. 

Key words: interacting particle systems, phase transition, metapopulation models, Allee 
' effect, mass migration, stochastic order, comparison with percolation. 
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1 Introduction 

<n : 

A metapopulation model refers to many small local populations connected via migrations 
, in fragmented environment. Each local population evolves without spatial structure, it 

| can increase or decrease, survive, get extinct or migrate from its site in different ways (see 

O ■ P2], for more about metapopulations). 

The most natural model for the evolution of a single population is the branching process, 
see [13] : birth and death rates depend on the number of individuals of the population, 
< ^ ' and the growth rate is density dependent. 

^ . If the birth rate is always larger than the death rate, if the population survives, it will 

increase indefinitely. If the birth rate is smaller or equal than the death rate, the population 
will get extinct almost surely (|27J). A more interesting situation is given by a birth 
rate larger than the death rate under a particular population size N, and smaller over 
that. The real environments observation suggests that this process is gradual, that is the 
growth rate decreases over a population size as population density increases. In some of 
our applications we suppose that over a fixed number of individuals N (the capacity of a 
site) the growth rate is null. 

Many biological phenomena may influence the dynamics of a metapopulation. 
Migration is one of the most important strategies that a species adopts to improve its 
probability of survival (see for instance [6], [15], [25]): when the population size is large 
one or more individuals leave the site where they are located to look for new resources in 
different sites. 

Other biological factors may favor the extinction of a species. One of them, the Allee 
effect, consists in an increase of the death rate when the density of individuals is small. 
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The reason is that at low density many factors (as difficulties in finding mates) cause a 
decrease of fecundity and an increase of mortality (see p], [9], [23], [25]). 
We simplify the real structure and we treat 4 metapopulation models from a mathematical 
point of view: we start from the easier one by adding a new biological phenomenon at 
each model. 

The mathematical models are interacting particle systems on Q = X^ d , where ICN: 
each particle represents one individual and on each site of Z rf there is a local population 
with capacity A (possibly A = oo), which evolves in different ways depending on the 
model. The local populations are connected via migrations of individuals, that is jumps 
of particles from a site to another one. 

In Section [2J we introduce the particle system, we give the main definitions and notations 
and we state the attractiveness results, crucial in the sequel for the existence of critical 
parameters and non-trivial invariant measures. Theorem 12. 11 the main result of a previous 
paper ((3J Theorem 2.4], inspired by [H]), gives necessary and sufficient conditions for at- 
tractiveness of a large class of particle systems. This simplifies many proofs, since in order 
to derive either if two processes are stochastically ordered or if a process is attractive we 
do not need to construct an explicit coupling for each model, but we only have to check 
inequalities involving the transition rates. 

In [20] and [22], the author considers a metapopulation model to investigate the roles 
of mass death (that is the death of all individuals in a local population) and spatial 
aggregation in the extinction of a species. In [2UJ he shows that, in presence of mass 
death, animals living in large flocks are more susceptible to extinction than animals living 
in small flocks: for that model, mass death can be an alternative to the Allee effect in 
raising to the exinction of a species. The new results in |22] involve the role of spatial 
aggregation, which may be either bad or good for survival in a model respectively with 
or without mass death. For these models the local population Allee effect was not taken 
into account. The model introduced, called non- catastrophic times model, is the following: 
for a fixed A < oo, on each site of T, d we may have up to A individuals: hence A is the 
capacity of sites. The transitions of the Markov process (rjt)t>o are 

■q t (x) r] t (x) + 1 at rate rj t {x)(p + %{ Vt ( y )=N} for <r] t (x) < N - 1, 

r]t(x) — > r]t(x) — 1 at rate 1 for 1 < r]t(x) < A, 

where y ~ x are neighbors. In other words each individual gives birth to another one on 
the same site with rate (p and dies with rate 1. An individual on site x gives birth to a new 
individual in a neighboring site with rate A/A only when the population at x has reached 
the maximal size A. There is a critical parameter for the capacity A of sites: 

Theorem 1.1 JUl Theorem 2] 

Assume that al > 2, A > and ip > 0. There is a critical value N c (\,ip) such that if 
A > N c {\,ip), then starting from any finite number of individuals, the population has a 
strictly positive probability of surviving. 

Starting from non-catastrophic times model, we propose 4 models to improve the un- 
derstanding of species dynamics. We want to investigate, for the first time in a model 
with spatial structure, the role of the Allee effect, the role of mass migration and their 
interactions. 
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In Section[3]we introduce Model I. This will represent our basic model with neither Allee 
effect nor mass migration. We begin with a system very similar to Schinazi's model: since 
a further step consists in adding migration of many individuals, we consider a migration of 
one individual to a neighboring site instead of a birth of a new individual. If 7V = 1 such 
a difference does not allow survival for the model with migrations, since no new births 
are possible, and the process gets extinct for any A: this is definitely not the case for the 
non catastrophic times model with N = 1, which is the contact process. If N is large this 
small difference does not change the behaviour of the model. 

This is the basic model and it must be as easy as possible (births and deaths on the same 
site and migrations from one site to another, all for at most one particle at time). For this 
reason we do not consider mass death, which is an additional complex factor. 
We take the birth rate larger than the death rate, but we fix a capacity N per site. A 
migration of one individual from a site x towards a nearest neighbor one is allowed only 
when the population on x reaches N. We prove that in some cases there is almost sure 
extinction and in others the species survives with positive probability: the key tool to 
prove survival is the comparison technique with a supercritical oriented percolation model 
(see [12]). 

In Section [5] we introduce Model II, that points out the key role of the Allee effect 
in species dynamics. Schinazi used mass death to prove that it can be considered an 
alternative to the Alle effect for extinction of a species. Since both the Allee effect and 
mass death improve the probability of extinction, in order to understand the role of one 
of them they should be considered separately. Here we want to show that a strong Allee 
effect (with neither mass death nor mass migration) is a key factor for the extinction. 
We add the Allee effect to Model I. Different probabilistic tools have already been used 
to illustrate the Allee effect like stochastic differential equations (see [10]). discrete-time 
Markov chains (see [2]) or diffusion processes (see [TT]). but none of these models has a 
spatial structure. 

In Model II each site has a capacity N, but the death rate is larger than the birth rate 
for small densities. Migration works exactly as in Model I. Theorem 14. II states that for all 
possible capacities, growth and migration rates, there exists an Allee effect large enough 
for the species to get extinct. It is proved through comparison with subcritical percolation. 

In Model III, introduced in Section we allow a migration of more than one individual 
at a time from one site to the neighboring one in a species affected by the Allee effect. We 
prove that mass migration might be the possible strategy of a species to reduce the Allee 
effect and improve its survival probability. 

When a local population size reaches N, a migration of a number of individual smaller 
than a fixed M is possible. In Model II, for an Allee effect large enough, the species gets 
extinct. In Model III, if N is large enough there exists M such that this is no longer true. 
A migration of large flocks avoids small densities in a new environment which are bad for 
survival. Indeed, by comparison arguments with oriented percolation, even if the Allee 
effect is the strongest one, if the species lives and migrates in flocks large enough, survival 
is possible (Theorem 15.1 p . 

In Section Owe generalize the previous models: in Model IV, instead of fixing a capacity 
N , we consider a slightly more realistic model. In all environments there is no maximal 
size, but a kind of self-mechanism of birth control such that the death rate is larger than 
the birth rate when there are more than iV individuals in a local population. A migration 
of one or more individuals is allowed from a site with more than N individuals towards 
a site with few individuals. We prove in Theorem 16.11 that in some cases we can have 
survival but on each site the population does not explode even if there is no capacity. 
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Namely on each site the expected value of the number of individuals is finite. In other 
cases the species gets extinct. 

Note that on each model instead of fixing the death rate equal to 1 and letting the birth 
rate vary (the most used approach), we consider the reverse but equivalent point of view 
in order to clarify our proofs, presented in Section [7J 



2 Background and tools 



The mathematical model is an interacting particle system (rjt)t>o on = X^ d , where 
X = {0, 1, . . . , N} C N and N denotes the common size (capacity) of the local popula- 
tions, if finite. The value r]t(x), x £ Z d is the number of individuals present in site x at 
time t > 0. We write O^r when we want to stress the dependency on the capacity N. 
When X is finite, which is the case of Models I, II and III, we refer to the construction 
in [T7J; when X is infinite, that is in Model IV, the state space is non-compact and a 
different construction is needed. The first examples of interacting particle systems with 
locally interacting components in non-compact state spaces have been introduced in |24| . 
One approach to construct these kinds of models has been developed in [TH], where the 
contruction was detailed for Coupled Random Walks but with small changes it can be 
generalized to many other processes. By using similar ideas, in [7J was stated a general 
existence theorem for reaction-diffusion processes, that we are going to apply in Model IV: 
in order to assure the existence of the process, some restrictions on the transition rates 
are required, as explained in Section El 

The process admits an invariant measure /j, if Pa(T]t G A) = ^(A) for each t > 0, A C Q, 
where P^ is the law of the process with initial distribution fx. An invariant measure is 
trivial if it is concentrated on an absorbing state when there exists one. The process is 
ergodic if there is a unique invariant measure to which the process converges starting from 
each initial distribution (see [HI Definition 1.9]). For any x,y £ Z d , we write y ~ x if y is 
one of the 2d nearest neighbors of site x. 

We introduce here a common infinitesimal generator C (we will be more precise on 
each model): it is given by 

= E E {piwtt&i) - f^+p-^ifis^vw) - m) 
+ E^ r S(^)(/(^)-/^)} 

where / is a local function, rj € f2, Sx,y ,k , S k and S~ k , where k > 0, are local operators 
performing the transformations whenever possible 

7](x) — k if z = x and n(x) — k E X, n(y) + k G X 
(Sx,y k v)(z) = { V(y) + k iiz = y and n(x) - k e X, V (y) + k e X (2.2) 
7](z) otherwise 

r)(y) + k if z = y and rj(y) + k € X 
rj(z) otherwise 



(S- k v )(z) = { V ^ ~ k if Z = V and V ^ ~ keX (24) 
y 1 7](z) otherwise 

P k , P.~ k are positive functions from X to R, and in our four models k = 0, 1 (particles are 
born and die one at a time). 
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We assume P Q = 0, that is the Dirac measure concentrated on the empty configuration 
Sq is a trivial invariant measure. The function 1^/ x , y \ represents the migration (jump) 
rate; a jump of more than one particle per time is possible. We call emigration from x 
a jump that reduces the number of particles on x and immigration a jump that increases it. 



There is a natural definition of partial order on the state space: 

V£,r/GQ, £ < no (Vs G S,£(x) < n(x)). (2.5) 

A process (rjt)t>o with generator C is stochastically larger than a process (£t)t>0 with 
generator £ if, given £o < % ; there exists an increasing Markovian coupling (£t,77t)t>o on 
state space fix!] such that 

P«°>"°)(£ t <r/ t ) = l 5 

for all t > 0, where p(&>> , ? t >) denotes the distribution of {S,t,Vt)t>o with initial state (£o>?7o)- 
In this case the process (Ct)t>Q is stochastically smaller than (rft)t>o and the pair (£ t , r]t)t>o 
is stochastically ordered, see [3] Section 2]. If C = C and there is stochastic order between 
two processes with ordered initial configurations, then the process is attractive, see [171 
Definition II.2.2]. 

Necessary and sufficient conditions for stochastic order and attractiveness in a general 
class of particle systems including the models defined by generator (I2.ip have been derived 
by [3 Theorem 2.4], which generalizes |144 Theorem 2.21]. Since ()2.ip involves neither 
births nor deaths depending on neighboring sites, this theorem can be restated as follows: 

Theorem 2.1 ([3, Theorem 2.4]) Given K G N, j := {ji}i<i<K, m := {mi}i<i<K , 
h := {hi}i<i<K, three non- decreasing K-uples in N, and a,/3,"/,6 in X such that a < 7, 
j3 < 5, we define 

K 



I a :=/f (j,m) = \J{keX:mi > k > 8 - p + j { } (2.6) 

i=i 

K 

I b := if (j, m) = \J{k G X : 7 - a + m > k > j,} (2.7) 

i=i 

K 

J c :=/f(h,m)= |J{A; G X : m 8 > k > 7- a + /i J (2.8) 
i=i 

K 

I d :=/f (h,m) = U{fcGl:<5-^ + mi > A; > (2.9) 



i=i 



An interacting particle system (r]t)t>o with transition rates {T* b , P^,P a k }{ a ,b,k&x} ^ s stochas- 
tically larger than a particle system (^)t>o with rates {T k . b ,P k ,P~ k }{ a bkeX} if and only 
if 

E p| + e f ^< E ^i + E r ^ ( 2 - 10 ) 

fceX:fc>5-/3+j! fceia leX:l>j! l&I b 

E ^+E f ^> E ^v' + E r 7,* 

/or aiZ choices of K , h, j, m, a < 7 and /3 < <5. 
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Remark 2.2 It is not possible that an infinite value for K , I a , If,, I c , 1^ results in the 
same rate inequality: therefore one restricts to take K smaller than the maximal change 
(birth, death or migration) of particles involved in a transition (see J3J Remark 2.5]). 

Remark 2.3 To prove Theorem \2.1\ following the approach of ]14l, we first show that 
Conditions i2.10\) - ft2.11\) are necessary. Then we construct a Markovian coupling which 
turns out to be increasing under \2.1U\) - [2.11\) (see [3, Section 3]). Hence if Conditions 
\2.10\) ~ l2.11\) are not satisfied it is not possible to find a coupling that preserves the order 
between the two processes. 

By taking two processes with the same transition rates, Theorem 1 2 . 1 1 states necessary and 
sufficient conditions for attractiveness. We use attractiveness of a process to construct a 
non-trivial invariant measure starting from an initial configuration rjo £ Q^r, where 

Tl N := {q £ n : r)(x) = N for all x G Z d }. (2.12) 



Remark 2.4 (J^ Proposition 2.7]) For processes with births, deaths and jumps of at most 
one particle per site, Conditions A2.10\) and \2.11\) reduce to 

F>} + f I tP <P} + if (3 = 5 and 7 > a, (2.13) 

Pj <P} if p = 5 and-/ = a, (2.14) 

P- 1 + fi )/3 >P^ + if 1 = aand5>p, (2.15) 

Pa 1 > p ^ if -i = a and 5 = (3. (2.16) 



Remark 2.5 By [3, Corollary 3.28], the sufficient condition still holds if we consider 
systems with more general transition rates 1^/ n (x, y) and P^ x ^(x), not translation 
invariant. In this case there is stochastic order if Conditions \2.19) - l2.11\) (resp. A2.13)) - 
$2.16}) if N = 1) are satisfied for each pair of sites (x, y) and configurations rj < £ with 
r]{x) = a, rj(y) = ]3, = 7, £{y) = S. 

Remark 12.51 will be used in some steps of the further proofs (for Theorems 13.21 and 14. ip , 
where in order to make a comparison with oriented percolation we will introduce systems 
with different transition rates in different space regions, so that they do not satisfy the 
hypothesis of Theorem 12.11 

Definition 2.6 For a process (r}t)t>o there is survival of the species if 

IP(M > 1 for allt>0)>0 (2.17) 

where \rjt\ denotes the number of individuals at time t and \rjo\ is finite. Otherwise the 
species gets extinct. If the process starts from an infinite 770 we say that the species gets 
extinct if the process converges to So . The convergence to 5q is intended that for any finite 
S C Z d , the probability that there exists to such that for all t > to, rj t (x) = for all x € S 
tends to 1. 
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3 Model I: the basic model 



We introduce Model I. We choose to fix a birth rate equal to 1 and to associate two 
parameters to death and migration rates. Given cf> and A positive real numbers, transitions 
are, for all x £ S, y £ 5, y ~ x (we follow notations in 112. IB ) 

r] t (x) ^r]t(x) + 1 
r]t(x) ^r]t(x) - 1 

ivt{x),vt{y)) -^(vt(x) - i,vt(y) + 1) 

The model has the following monotonicity properties 

Proposition 3.1 Let (£t)t>o, (vt)t>o be two processes with respective parameters (0i, XN) 
and (4>2,X,N) such that (pi < <p2- Then (£,t)t>o is stochastically larger than {r]t)t>o, and 
(Vt)t>o is an attractive process. 

The key for attractiveness, which is a consequence of the stochastic ordering when (f>i = (p2, 
is that there are births, deaths and migrations of at most one particle per time and the 
migration rate from r)t(x) to rjt(y) is non-decreasing in rjt(x) and non-increasing in rjt(y). 

Corollary 3.2 Given (f]f)t>o such that t/q = £, then 

P(|»7t I > 1 for all t > 0) 
is non-increasing in (j) for each ^£(1. 

Remark 3.3 There is no stochastic order between systems with different values of N or 
A. Indeed in these cases conditions of Theorem \2.1\ are not satisfied. 

The first result corresponds to Theorem 11.11 for non catastrophic times model and it is 
proved in a similar way. 

Theorem 3.1 Suppose d > 2, A > and < 1. There exists a critical value N C (X,4>) 
such that if N > N c (X,(p), then starting from r]Q € such that 1 770 1 > 1 the process has a 
positive probability of survival. Moreover if rjQ G £1^ the process converges to a non-trivial 
invariant measure with positive probability. 

Proof . We skip it, since the result is a corollary of Theorem 15.11 We can get an easier 
proof that the process has a positive probability of surviving by slightly modifying [221 
Proof of Theorem 2] . The differences are that we consider a migration instead of a birth 
from x to y ~ x and the migration rate from x to y is non increasing in /?t(y). Such 
changes are not relevant for the proof. □ 

As we can expect aggregation is good for Model I, as in non catastrophic times model. 

Remark 3.4 If N = 1 the process dies out, since each individual can only migrate or die. 
This suggests that an increase of N is good for the survival of the species. However, by 
Remark \3.3\ there is no monotonicity property with respect to N. 

If we fix the capacity N, we prove that there is a phase transition also with respect to the 
death rate </>. 



at rate P m{x) = n t (x)t {rit{x)<N} , 
at rate P~K = firjt (x), 



(3.i; 
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Theorem 3.2 For all A > 0, 1 < N < oo, there exists c (A,A) < 1 such that, if 
4> < 4> C (X,N) the process starting from rjo with 1 < [770 1 < 00 has a positive probability of 
survival and if (ft > cj> c (\,N) the process dies out. Moreover for rjo E J7jy if 4» < </>c(A, A) 
the process converges to a non-trivial invariant measure with positive probability. 

We prove it in three steps, in Section Yl. 1.21 First (Step (i)) we find </>,!, (A, A) small enough 
to have survival: by Proposition 13.11 the process survives for each <f> smaller than <f>\{\, A). 
Then (Step (ii)) we prove that the process dies out for all A, A by taking <f> > 1 if it starts 
from a finite initial configuration and by taking (j) > 1 if it starts from 7/0 £ Finally 
in Step (iii) we use Corollary 13.21 to obtain the existence of a critical parameter C (A, A). 
Figure [U sketches the phase diagram in the (A, <^>) plane. 



Extinction 



1 



c (A,JV) 



Survival 



Figure 1: Phase diagram of Model I for a fixed N > 1 and finite initial configuration: by Theorem 
I3.2l there exists a critical curve </> c (A, N) which converges to as A goes to zero; it is always smaller 
than 1. We conjecture that <f> c (\, N) is monotone and as A goes to infinity converges to a value 4>m 
depending on the capacity N of the model which is strictly smaller than 1 for each N < 00. 

The model admits a phase transition with respect to the death rate <fi for each N > 2, 
while the same process without migrations dies out almost surely. The effect of a migration 
is to move an individual from a site in state N, where there is no possibility to give birth, to 
a site with less than N individuals, where it may reproduce itself. Therefore even if there 
is no monotonicity with respect to A (cf. Remark I3.3H , this suggests that an increase of A 
is good for survival. Contact interactions and migrations work in a similar way, but small 
differences are present. From a mathematical point of view an increase of the migration 
rate does not favor ergodicity. 

4 Model II: the Allee effect 

We translate the Allee effect into mathematical terms for a metapopulation model. As in 
Model I, we fix a capacity N for all sites, but we assume the death rate larger than (or 
equal to) the birth rate when the density is small. Namely, fix a positive integer Na < A 
and positive real numbers (j), A and 4>a > 1; the transitions are, for all x € S, y £ S, x ~ y, 
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referring to notations in (|2.1|) 



(Vt{x),r] t (y)) -+(m( x ) - 1 ,Vt(y) + 1) at rate — r. 



r] t (x) -^r]t(x) + 1 
rj t (x) -+r]t{x) - 1 



at rate P 



at rate P 




(4.1) 



We assume a > 1 and a > <f>) in other words if r] t (x) < Na then the death rate (pAVt( x ) 
is larger than (or equal to) the birth rate rj t (x) because of the Allee effect. If rjt(x) > Na, 
the most interesting situation is given by a death rate <f>rjt(x) smaller than or equal to the 
birth rate rj t (x), that is < 1. If either > 1 and r/o is finite or > 1 and tjq £ &n the 
species gets extinct as proved in Theorem 13.21 If Na = (no Allee effect) or Na = N 
(death rate always larger than birth rate) there is only one death rate and we are back to 



Since only births, deaths and migrations of at most one particle are allowed and the 
migration rate from r)t(x) to rjt{y) is non-decreasing in r]t(x) and non-increasing in r/t(y), 
attractiveness conditions are satisfied. One proves in a similar way that Proposition 13.11 
still holds for Model II either with respect to a or 0. Namely 

Proposition 4.1 Let (£t)t>o an d (i}t)t>o be two Model II type processes with respective 
parameters (0i, 0,4,1, A, N, Na) and (02, 4>a,2, A, N, Na) such that 0i < <p2 and 4>a,i < 4>A,2- 
Then (£,t)t>o is stochastically larger than (rjt)t>o and (rjt)t>o is attractive. 

The corresponding Corollary 13.21 holds in a similar way. 

We prove that the Allee effect changes the behaviour of the system: for any possible 
capacity and migration rates there exists an Allee effect large enough for the species to 
get extinct. 

Theorem 4.1 Assume (f>A > 1 and let 4> C (X,N) be the critical parameter introduced in 
Theorem\3M Then for all X > 0, < N < oo, < N A < N: 

i) if (f> < <p c (\,N), there exists $5(0, A, N, Na) such that if 4>a > 0^(0, A, N, Na) the 
species gets extinct for any initial configuration r/o £ and if 4>a < 4>^~{4>, A, N, Na) 
the species has a positive probability of survival; 

ii) if C (A, N) < (< 0a) ; the species gets extinct for any initial configuration tjq £ fijv- 

This corresponds to the biological idea that random fluctuations, which are present on 
each local population, plus the Allee effect doom even a very large population. 
The phase diagram of Model II depends on 0a- Proposition 14. II is not enough to construct 
a detailed phase diagram, but it gives some information in this direction. Since for any 
and A there exists 4>a large enough for the species to get extinct, one can choose 0a large 
enough to reduce the survival region in the (A, 0) plane of Figure [T] for such fixed 0a- 

In order to model the Allee effect, we require 0a > 1 and < 1. Note that from a 
biological point of view we just need 0a > 0, but if either 0a > > 1 or 1 > 0a > 0, by 
monotonicity arguments we can work as in Model I. 

From a mathematical point of view, it would be interesting to investigate a model where 
and 0a play symmetric roles, that is 0a < 1 and > 1. For fixed N, Na and A we 
prove that there is no 0a such that there is survival for all and no such that there is 
extinction for all 0a- 



Model I. 
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Theorem 4.2 For all 1 < N A < N , A > 0, 

i) for each <f> > 1 i/iere exists (f)^(X, Na, N,<j)) such that, if 4>a < ^(X, -^4, N, <j>), 
the process survives for any initial configuration 770 such that \rjo\ > 1 with positive 
probability; 

ii) for each <f>A < 1 i/iere exists <p c (X, Na, N, 4>a) such that, if 4> > <p c {\ N A , N, 4>a), the 
process dies out for any initial configuration t]q G Q^. 

5 Model III: Mass migration as Allee effect solution 

We have already observed in Model I that a migration of a single individual is good in 
absence of the Allee effect. The model without migrations dies out, but if we add a pos- 
sible migration of one individual there is a positive probability of survival. In Model II, 
anyhow, a single individual migration may not be enough: even in the supercritical region 
of 4> in Model I there exists an Allee effect strong enough for the species to get extinct. 
Which strategy may a species adopt to reduce the Allee effect? 

We show that, at least in theory, migrations of large flocks of individuals improve the 
probability of survival for any Allee effect. A migration of many individuals in a new envi- 
ronment improves the probability of a successful colonization avoiding a small density in 
that new environment which is influenced by the Allee effect. 

We introduce positive parameters 4>a, <fi, N A , N such that < Na < N, <f) A > 1, <j> > 
and we take birth and death transitions as in Model II, but more general migration rates: 
given M € N, < M < N, y ~ x the transitions are 

7fc(x) -^Vt(x) + 1 at rate P^ x) = Vt(x)l{ Vt ( x )< N -i}, 

n t (x) -*th(x) - 1 at rate P~^ x) = Vt(x)((t>Ai{ Vt ( x )<N A } 

+ 4>^{N A < Vt (x)}), (5-1) 

{Vt(x), r) t (y)) ^(rit(x) - k, r/ t {y) + k) at rate ^^ t ( x)My ) = ^{vt(x)-k>N-MMy)+k<N} 

for 1 < k < M. In other words if k G {1, 2, ... , M} individuals try to migrate from x to y 
but if rjt{y) + k > N, the migration does not happen. Notice that if rft(x) < N — M the 
migration rate is null: individuals try to migrate only when there are more than N — M 
individuals on a site. From a biological point of view, this means that when there are few 
individuals resources are enough for all and there are no reasons to migrate. When rj t (x) > 
N — M there is a positive probability of migration and the number of individuals that may 
migrate is increasing with the population size. If rjt(x) = N — M + 1 we allow a migration 
of at most 1 individual from x to a nearest neighbor site, when rj t (x) = N — M + 2 we allow 
a migration of either 1 or 2 individuals with rate A and so on. If n t {x) = N = (N — M) + M 
we allow a migration of 1, 2, . . . to the largest flock of M individuals, where each migration 
occurs with rate A. 

First of all we notice monotonicity properties. 

Proposition 5.1 Let (£t)t>o an d (??t)t>o be two Model III type processes with respective 
parameters 0^,1, X,N, Na) and (cp2, 4>A,2, A, N, Na) such that 4>i < 4>2 and 4>a,i < 4>A,2- 
Then (£t)t>o is stochastically larger than (rjt)t>o, o,nd (rjt)t>o is attractive. 
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The corresponding Corollary 13,21 holds in a similar way. 

In Model II we showed that a strong Allee effect dooms even a very large population 
with a large migration rate. The strategy that the species may adopt to reduce the Allee 
effect is to increase the number of individuals which migrate: we prove that we can take 
a population size N and a maximal migration flock size M large enough for the species to 
survive for any Allee effect. 

Theorem 5.1 Let d>2. For all X > 0, N A > 0, 

i) if (ft < 1 there exists N c ((ft, A, Na) such that for each N > N c ((ft, A, Na), there exists 
M(Na) so that the process starting from t]q with \rjo\ > 1 has a positive probability of 
survival for each (ft a < oo. Moreover ifrjo G fijv the process converges to a non-trivial 
invariant measure for each eft a < oo; 

H) if (ft > lj the process gets extinct for all N, X, (ft a > lj M and for any finite initial 
configuration. If r/o £ Qjv is not finite the process gets extinct if (ft > 1. 

Remark 5.2 T/ie proof of i) (see Section^ states that in order to have survival we can 
take M(N A ) = N A + 1. If N A = 0, this gives M(N A ) = N A + I = I, only a migration 
of one individual is possible and the process reduces to a Model I type process: therefore 
Theorem \3.1\ is a particular case of Theorem \5.1[ 

Notice that N c ((ft, X, Na) does not depend on (ft a- This means that even if the Allee effect 
is the strongest one, if the species lives and migrates in flocks large enough, survival is 
possible. 

Since there are many parameters the phase diagram is not easy to construct; nevertheless 
Proposition 15.11 suggests that one can choose N and M large enough to extend the survival 
region in the (A, </>)-plane for fixed cf>A, N and M. 



6 Model IV: ecological equilibrium 

Real natural environments do not have any a priori bound on the population size, but there 
is a kind of self-regulating mechanism that does not allow an "explosion" of the number 
of individuals per site. Ecological equilibrium has been introduced in [3] for restrained 
branching random walks (on a connected, non-oriented graph X with bounded geometry) 
with transition rates: 

7](x) — > 7](x) + 1 at rate r](y)p(y, x)c(j](x)) 

y 

rj(x) —¥ rj(x) — 1 at rate rj(x) 

where c : N — > IR + is a non-increasing function and P = (p(x,y)) Xy y^x is a stochastic 
matrix such that p(x, y) > only if x ~ y. The idea is that some restrictions on branching 
random walks birth rates, given by the non-increasing function c(-) of the number of indi- 
viduals, provide survival within nonexploding populations. In particular, one interesting 
consequence of [H Proposition 1.1] is that one can find a function c such that the process 
survives but limsup^^ W°{rjt{x)) < oo uniformly for any bounded r)Q £ O and x £ X. 
We show that a similar mechanism leads to a similar conclusion on different systems. In- 
stead of taking births on neighboring sites as in [3], we consider a non-increasing birth 
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rate in the same local population but we add migrations when the number of individuals is 
larger than a fixed value N. This means that the restriction on birth rate does not change 
the migration rate: this is not the case for the restrained branching random walk, where 
births in a new site (which play the same role as migrations in Model IV) depend on the 
local population density. 

We suppose that in our environment there is no maximal population size as in previous 
models, and the birth rate is always positive. We also assume that, when the population 
size is larger than N, the death rate increases faster than the birth rate, hence the growth 
rate is negative. 

In order to simplify notation and proofs, we work on a modification of Model I. Namely, 
given positive real values (ft, (ft, we take the following transitions, for each x G S, y G S, 
x ~ y, 

Vt(x) ->ri t (x) + 1 at rate P^ t{x) = n t {x), 

Vt(x) -HftO) - 1 at rate P~ l {x) = Vt(x)(cftt {r]t(x) < N} + (ftl {N<Vt(x)} ) , 

(6.1) 

1 A 

{Vt(x),Vt(y)) ->(vt(x) - l,Vt(y) + 1) at rate —T^^ = ^■{ m (x)>N, m ( y )<N}- 

This means that when the population size rjt(x) is larger than TV the death rate (ftrjf(x) 
is larger than the birth rate rjt(x). A migration is allowed from a site with more than TV 
individuals to a site with less than TV individuals. Since we are working without any a 
priori bound, we refer to construction techniques in non compact cases and we restrict the 
state space to Q, C fi (See [3 Chapter 13]), where 

Q := {77 G $7 : f]{x)a(x) < 00} 

and (oc(x)) x£ %d is a positive sequence such that X^ez d a ( x ) < 00. Sufficient conditions 
for existence and uniqueness of the process given in [7J Chapter 13] are satisfied: 

Lemma 6.1 There exists a unique Markov process with state space 0, generator \2.1\) 
and rates \6. 

Since births, deaths and migrations involve only one particle and the migration rate is 
non-decreasing in r]t(x) and non-increasing in r]t(y) the process is attractive as in Model 
I and a monotonicity property (see Proposition I3.2|) holds in (ft and in eft for each initial 
configuration r]Q £ fl. We prove that in some cases the process survives but does not 
explode, that is it does not die out, and the expected value on each site is finite. 

Theorem 6.1 Let tjq G 0, n for some n G N (so that 770 G Q,). For all A > 0, (ft > 1, 

i) for each 1 < N < 00 there exists (ft c (X,N,(ft) > such that if (ft < (ft c (X, N,(ft) the 
process has a positive probability of survival and if (ft > </> c (A, N, (ft) the process dies 
out; 

ii) for each (ft < 1 there exists N c (\,<ft,(ft) > such that if N > N c (X,(ft,(ft) the process 
has a positive probability of survival. 

If the process survives, there exists C n < 00 so that lim E(r] t (x)) < C n for each 

t— >oo 

x G Z d . 
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Note that the constant C n depends on the initial configuration. Since the migration rate 
does not depend on the local population density we are not able to find such a constant C 
independent of the initial configuration, which was the case for the model treated in [3]. 



Remark 6.2 In a similar way one can consider a Model III type process without any 
a priori bound by adding a death rate 4>7]t(x) when the number of individuals in a local 
population is larger than N . By comparison arguments, even if a strong Alice effect is 
present a mass migration of large flocks of individuals leads to the survival of the species 
but the local populations do not explode. 



7 Proofs 

We first recall a classical result involving random walks on a finite interval. Let r\, r2 € N 
and (Xt)t>o be a discrete time random walk on {r±,ri + 1, . . . ,r% = r\ + n} such that 

i — > i + 1 with probability p, i E {r\, . . . ,r% — 1}, 
i — > i — 1 with probability q, i £ {ri + 1, . . . , 7*2}. 

We interpret this random walk as a game which ends when Xt reaches either r% or 7*2 , that 
we call respectively the ruin of the first and the second players. 

Lemma 7.1 (Ruin Problem Formula, \19V (4-4); Section 1-4])- Let P r2 (j) (resp. P Tl {j)) 
be the probability that the random walk starting at j £ {r% + 1, . . . , 7*2 — 1} reaches state r2 
before state r\ (resp. state r\ before r2). Then 

1 _ (q/nV- ri 



7.1 Model I 

7.1.1 Proof of Proposition I3TT1 

We prove that if 770 < £o, then rjt < £t for each t > a.s. This is an application of Theorem 
12.11 since there is a change of at most one particle per time, we check conditions in Remark 
12.41 The transition rates are given by (13. with 4> = <f>2 for the process (rjt)t>o and 4> = 4*1 
for (£t)i>o- Conditions (|2.13p and (|2.14p are: given 77 < £, if r](y) = and 77(2) < £(2;) 

r l{y)^{r ! (y)<N-l} + ^{r](x)=N, V (y)<N} <C(y)^{$(y)<N -1} + ^{Z(x)=N£(y)<N} , 
r](y)^{ V (y)<N-l} <£(y)%(y)<7V-l}- 

Since 77 < £ and 77(2) = N imply £(x) = N and since l{ v (x)=NMy)<N} < %(aO=Ar,£(2/)<iV} 

if r/(y) = £(t/), the conditions are satisfied. 

Conditions (|2.15|) and (j2. 16[) are: if rj(x) = and 77(7/) < 

(p2V( x )^{v(x)<N-l} + ^{ v (x)=N, v (y)<N} ><Plt,{ x )^{t(x)<N-l} + ^{£(x)=N£(y)<N} , 

which hold since 4> 2 > 4>i and t {v{x)=N ^ y)<N} > %( ffi )=Ar,5( 2 /)<iV} , because rj(x) = £(x).D 
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7.1.2 Proof of Theorem ET21 



We prove it in three steps. In Step (i) we find cft],(X,N) small enough to have survival; in 
Step (ii) we prove that the process dies out for all A, iV by taking > 1 if it starts from a 
finite initial configuration and by taking (ft > 1 it it starts from rjo G J7jy, and in Step (Hi) 
we get the existence of a critical parameter by monotonicity. 

i) We follow the idea in [21] by using the comparison technique with oriented perco- 
lation (introduced in [5]) explained in [12]. Here and in the subsequent proofs we think 
of the process as being generated by the graphical representation, see [T2] for such a con- 
struction. Suppose d = 2. The proof in higher dimension is similar but the notation is 
more complicated. Denote by 

ei = (1, 0), M = {(m, n) G Z 2 : m + n is even} 

B = (-4L,4L) 2 x [0,T], -B m ,n = {2mLe\,nT) + B (7.1) 
/ = [-L, L] 2 , I m = 2mLei + / 

where L and T are integers to be chosen later. In other words B m ^ n is the cube that we 
get by applying a translation of (2mLei,nT) to B and I m the square we get by applying 
a translation of 2mLei to /. Roughly speaking, the idea consists in constructing boxes 
large enough so that with large probability the species survives inside a box, and then to 
compare this evolution with an oriented percolation model. 

Let (r] t )t>o be the process defined by generator (12. 1|) with rates (|3.1f) . We consider a 
modification rff l ' n of rj t : the process (^™'™)t>o is constructed through the graphical repre- 
sentation of r\t in B mt n, but rj t (x) = for all x ^ B m ^ n and t > 0. Let m Xj y with y ~ x be 
the Poisson process with rate \/(2d) associated to a migration from x to y. A migration 
from x belongs to the graphical construction in B m ^ n if x G B m)Tl : therefore an immigra- 
tion to B mtn from a site y ^ B m)Tl cannot happen for r/™' n , but we still consider the arrows 
of emigrations from B m ^ n . Their effect is the death of one individual on the boundary of 
B m ^ n . If t]q(x) = rj^ ,n [x) = ]l{ y }(x) for some y G B mjn , rj t > r(^ ,n by Remark 12.41 since if 
x ^ B mtn then ?7™' n (x) = 0, otherwise conditions in Remark 12.41 are satisfied for each pair 
of sites (x,y) (see also Remark l2.5p . 

We say that (m, n) is wet if rj^ 1 ' 71 starting at time nT with at least one individual in I m 
is such that there is at least one individual in I m -i and one individual in I m +\ at time 
(n + 1)T. Otherwise the site is dry. The event G m ^ n := {(m,n) is wet} is measurable 
with respect to the graphical construction in B m ^ n : we prove that we can choose L and 
T such that the probability of a site (m, n) to be wet can be made arbitrarily close to 1 
if <f> is small enough. By translation invariance it is enough to show it for (0,0). We call 
Vt' := (t, we fix L > and we prove that for each e > there exists T and <fi such that 

P((0,0) is wet) > 1 - e, (7.2) 

that is that if there exists one individual in a site G Io =: / there is at least one 

individual both in I\ and i_i with large probability. 

In order to prove it for (j> small enough, we begin by showing that it holds for a process 
with (ft = inside B: let P(£i G •) denote the law of such a process. This means that each 
individual in box B survives forever. 

We choose a preferential path (i, j), [i + 1, j), . . . , (L,j), (L + 1, j): we prove that there 
exists T large enough so that the abscissas of the rightmost and leftmost particles are 
respectively larger than L and smaller than —L with probability larger than 1 — e, since 
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this is one possibility for the site (0, 0) to be wet. 

We wait until in we have a stack of N individuals and then wait for the first migration 
to (i + 1, j) (we may have to wait to pile up again after a migration in a wrong direction). 
After the first migration from to + we wait for another migration from (i + l,j) 
to (i + 2, j) and so on so that in a finite time we reach (L + 1, j). A similar idea works for 
the leftmost particle. We conclude that if 4> = for all e > 0, A > 0, 1 < N < oo there 
exists T = LT such that 

p((0,0) is wet) > 1 - e/2. (7.3) 

Now we prove (17. 2p for (j> small enough. Let Al = Al(<P, N) be the time of the first death 
on the finite box (— 4L,4L) 2 . If Al > T := LT, for each e > we can take (j) > small 
enough for 

(0,0) is wet) >P((0,0) is wet\A L > T^F(A L > T) > pY(0,0) is wet) e ~^ N ^ L f T 
>1 - e. 

Hence for all e > 0, L > 0, A > 0, 1 < A < oo there exists T and <t>l(X, A) > such that 
if <t> < (j>l(\, N) then fl£2) holds. 

By comparing the process with an oriented percolation process, the existence of an infinite 
path of wet sites corresponds to the existence of individuals at all times and for e small 
enough percolation occurs (see [E]). By monotonicity (Proposition [3J]), the process sur- 
vives for any (j) < </>c(A, A). 

ii) Let £t be a continuous-time Galton- Watson process without spatial structure starting 
from |?7o| < Co individuals. We couple the total number of particles of the two processes. 
Each individual in both processes breeds at rate 1 (except for r\t when the full carrying 
capacity of the site is reached) and dies at rate 4>. Since we are interested in the total 
number of particles, migrations do not count in this coupling. Therefore \r]t\ < £t for all 
t > 0. If Co is finite and <p > 1 then the Galton- Watson process gets extinct; this implies 
that r\t dies out for any <j) > 1. 

Assume now that 770 £ ^jvj we prove that the process gets extinct when cf> > 1. By 
translation invariance, for each t > 

^E(?y t (x)) =K(Cr] t (x)) = K(r] t (x)l {vt(x) < N ^ 1} - <j)r] t (x) + ^ l {vt(y)=NMx)<N} X/(2d) 
~ 1 {vt(x)=NMy)<N} x /( 2d )) = E (^(a ; ) ]1 {»;t(a : )<JV-i} - <t>m(x)) 

y~x 

<(1 - 0)EMx)). 

and by Gronwall's Lemma the process converges to uniformly with respect to x. By 
Corollary 13.21 the process dies out for each (j> > 1. 

Hi) The claim follows by Steps i), ii), and Corollary 13.21 Starting from t]q € fJjv, the 
existence of the upper invariant measure follows from attractiveness, and it is non-trivial 
by Step i). □ 



15 



7.2 Model II 



7.2.1 Proof of Theorem H/Ll 

ii) Since 4>a > (f>, Model / is stochastically larger than Model 11. If <j) > </> c (A, N), both of 
them die out by Theorem 13. 21 

i) Assume <j) < 4> C (X,N) (< 1 by Theorem I3~2"j) . We follow the idea in [2^1 Theorem 4.4] 
and we compare the system with a subcritical percolation process. We prove i) when d = 2 
in order to simplify the notation (the same proof works for all d > 1). Let (rjt)t>o be a 
process with generator (12 .If) , rates (14. ip and r/o £ &n- We define 

' A = [-2L,2L] 2 x [0,2T]; £> = [— L, L] 2 x [T, 2T] 

( C b = {(x,y,t) eA:t = 0} 

C s = {(x,y,t) £ A : \x\ = 2L or \y\ = 2L} { > 

k C = C b UC s = {(x,y,t) £ A : \x\ = 2L or \y\ = 2L or t = 0}. 

where T is a time to be fixed later. 

In other words C is part of the boundary of the space-time region A, which contains the 
smaller region B. We construct a percolation process on M = 1? x Z + starting from 
(Vt)t>o- We consider for each (m,n,k) £ AT & modification r ]^ n,n,k of rj t - the process 
(j]™' n ' )t>o is constructed through the graphical representation of r/ t in A + (mL, nL, kT), 
but rj™' n ' k (x) = N for all x G {mL,nL) + (-2L,2L) 2 , t < kT and x (£ (mL,nL) + 
(-2L,2L) 2 for all i > 0. Therefore an emigration from A + (mL, nL, kT) cannot happen 
and an immigration from a site y on the boundary of (mL,nL) + [— 2L,2L] 2 after kT is 
always possible with rate A. By Remarks 12.41 and 12. 5( rjt < f?™'™' for all m, n, k and t > 0, 
since if x ^ (mL, nL) + (— 2L, 2L) 2 then ry™' n (x) = JV, otherwise conditions in Remark 12.41 
are satisfied for each pair of sites (x,y). 

We say that a site (m, n,k) £ J\f is wet if there are no individuals for the process r]™' n,k in 
B + (mL, nL, kT). A site is dry if it is not wet. 

We show, through a series of lemmas, that the probability of a site to be wet is as large as 
we want by taking 4>a large. By translation invariance we prove it for (0, 0, 0) and we denote 
^o,o,o ._ ^ q ^ ^ ^ First of all we prove that there exists a time S at which 
with large probability there is at most 1 individual per site on (— 2L,2L) 2 (Lemma 17. 3p . 
After S, there exists a time T such that there are no individuals in (— 2L, 2L) 2 with large 
probability (Lemma I7.4p . Therefore with large probability the only possibility of having 
one individual in B is that an emigration from the boundary after time T reaches [— L, L] 2 
before 2T: the last step consists in proving that such an event has small probability. 
We first introduce an auxiliary process whose transitions are not translation invariant: 

Lemma 7.2 Let (£t)t>o be a process with only birth and death rates: 

P)(x) = l { i< N -i } (l + X); P7\x) = l{<t)A\i<N A } + ^{N A <i}) ifx £ (-2L,2L) 2 

(7.5) 

and £t(x) = N for all x ^ (— 2L,2L) 2 , t > 0. Then (£, t )t>o is stochastically larger than 
(£,t)t>o- 

Proof. Both £t an d £f are equal to for each t > outside (— 2L,2L) 2 . By Remark 
\2.5\ we check Conditions in Remark 12.41 for each pair of sites (x, y) with either x or y in 
(— 2L,2L) 2 . If x £ (— 2L,2L) 2 , (Ct( x ))t>o is a birth and death process whose birth rate is 
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the original one plus the largest immigration rate on £t(x), and whose death rate is the 
original one plus the smallest emigration rate on £t(x), which is null. For each r\ € f2 



P v(l) + T l(x)Mv) ^M^NA^AVix) + l {NA<rt{x)} (jyq{x) = P~ {x) . 
then all conditions are satisfied. □ 

Lemma 7.3 For all e > 0, L there exists S > and 4>a such that 

P{G L (S)) > l-e/6 (7.6) 
where G L (S) = {£ s {x) < 1 for each x G (-2L,2L) 2 }. 

Proof . We prove (|7.6p for (^)t>o with law P(£ t € •). By monotonicity (Lemma IT.2[) it 
will be true for (£t)t>o- For all e > and L we take S large enough so that the number of 
visits H x to of £ t (x) before S satisfies 

p <"- = °» £ WW 

If there is at least one visit, we consider 

K oo 

^p(? 5 (x) > l|#f = k)f(H^ = k)+ £ ¥ {Zs(*) > = = fc) (7-8) 

fc=l k=K+l 

By taking if large enough the second sum (in which there are more than K hits to 0) is 
as small as we want. There are at least two individuals in a site after the i-th visit to 
only if the exponential clock Bi ~ Exp(l + A) (birth rate if (, t (x) = 1) rings before the one 
of D{ ~ Exp{<pA) (death rate if £t(x) = 1). Therefore for all e > 0, L and K we can take 
4>a large enough for the first sum in (|7.8p to be smaller than 

gF( 3i e {1, 2, . . . *} : B, < A) < K 2 ^^ < ^ (™) 

By (17.7P and (I7.9P for all e > 0, L there exists S and cpA large enough for 

P((G L (,S)) C ) < (4L) 2 sup P(| s (x) > 1) < e/6 (7.10) 

xG(~2L,2L) 2 

and the claim follows. □ 



Lemma 7.4 For a// L ; e > i/iere exists S and <f>A such that 

P(G L (S + S)) > l-e/3 
where Gl(S + S) = {£, s ,g(x) = /or eac/i x £ (— 2L,2L) 2 } and S is given by Lemma 7.3 
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Proof . If Gl(S) holds, we take S small so that there are neither births nor immigrations 
from the boundary C s between S and S and (pA large so that all individuals in (— 2L, 2L) 2 
die before S with large probability. Namely, given D ~ Exp((/)a) and B ~ Exp(jl + 

A)(4L — l) 2 ^, for all e > 0, L there exists 5 small and (Pa(S) large enough for 

W(G L (S + S)\G L {S)) >F(D < S)( 4L -V 2 W(S < B) 

>(1 - exp(-^)) (4L " 1)2 exp(-(l + A)(4L - ifS) 
>1 - e/6. 

If T = 5 + S 1 , given by the two previous lemmas, 

F((G L (T)) C ) <e/6 + e/6 = e/3 (7.11) 

and the claim follows. □ 

Therefore Ct(x) = for each x G (— 2L, 2L) 2 with large probability. Since Po 1 = 0, the only 
way to get an individual in [-L, L] 2 between times T = S + S, given by the two previous 
lemmas, and 2T is that a migration from y £ Ct = {y = (2/1,2/2) £ -4 : 1 2/1 1 = or 1 2/2 1 = 
2L} gives birth to a chain of individuals which reaches [-L, L] 2 in a time smaller than T. 
Suppose that £t(y) = iV for all y £ Ct and i G [T, 2T]. By monotonicity it will be true for 
any smaller configuration. We fix K large so that the number of emigrations -Er,c T horn 
Ct to (— 2L, 2L) 2 from time T to 2T is larger than K with probability smaller than e/3. 
After one migration, with probability smaller than (1 + X)/{<j)A + 1 + A) there is a new 
birth or a new immigration at x before the death of the individual. If the number of such 
migrations is smaller than K, by taking <pA large enough 

P((0,0,0) is dry \G L (T), E T)Ct <K)< ^\\+ x < (7.12) 

By (ITTTTD and (17121 we get 

P((0,0,0) is dry) <P((0,0,0) is dry \G L {T)) + e/3 

=P((0,0,0) is dry \G L (T),E T ,c T > K)P(E T ,c T > K) 

+ P((0,0,0) is dry \G L (T),E T ,c T < K)F(E T ,c T < K) + e/3 
<e/3 + e/3 + e/3 = e. (7.13) 

Now we construct a dependent percolation model such that the probability of a site to 
be wet is as large as we want. For all (m,n,k) and (x,y, z) in Af such that k < z and 
the intersection between (mi, nL, kT) + A and (xL, yL, zT) + A is not empty we draw an 
oriented edge. Notice that the probability of a site (m, n, k) to be wet depends only on 
the existence of a path of individuals within (mL, nL, kT) + A; since each block intersects 
only a finite number of other blocks, there exists K such that all sets of sites in M with 
distance larger than K are independently wet. Here the distance is the minimal number of 
edges (without orientation) connecting two sites. Therefore this is a dependent percolation 
model with finite range of interactions. 

By monotonicity, the probability of having an individual in metapopulation model (rjt)t>o 
in (mL, nL, kT) + B is smaller than the probability of the existence of a path of dry sites in 
the percolation model with endpoint (m,n,k) starting from (y, z, 0) for some (y,z) 6 Z 2 . 
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By working as in [26} Proof of Theorem 4.4], for any given site x G S there exists a random 
time T x a.s. finite after which there will never be any individual. Let A be a finite subset 
of 7L d and Ta '■= max{T x ,x G A}. By monotonicity, Ta may be chosen uniformly in the 
initial configuration 770. Given 770 € at, let V be the invariant measure lim 5 m T(t) (where 

t— >oo 

T(t) is the semi-group of the process), which exists by attractiveness. For each finite set 
A C Z d 

V(£ G SI : > for some x G A) = 0. 

Since 17 gives null probability to each set of configurations with at least one individual, it 
concentrates on the empty configuration, that is V ~ Jo, and ergodicity follows. 

7.2.2 Proof of Theorem 14^21 

i) We work as in proof of Theorem 13.21 with the same notations: we suppose d = 2, we 
use (|7,ip in order to make a comparison with an oriented percolation model and we define 
for each (m, n) a modification r]™' n of the process in the same way. A site (m, n) G Z 2 is 
wet if r] t ' starting at time nT with at least one individual in I m is such that there is at 
least one individual in l m ~\ and one individual in I m +i at time (n + 1)T. By translation 
invariance we work on ^ := 7y t '°. We will prove the analogue of (|7.2p . 

We start with one individual at x = G I and we choose a preferential path (i + 
. . . , (L,j), {L + 1, j): if there exists T such that the abscissas of the rightmost and 
leftmost particles of £t are respectively larger than L and smaller than —L at T then site 
(0,0) is wet. We begin by working with 4>a = and call P(£f G •) the law of the process 
in this case. 

We fix L > 0. We wait until in we have a stack of Na individuals: since A = 

{Na, Na + 1, • • • , N} is an absorbing set (because 4>a = 0) after a finite time the local 
population size reaches iV and migrates to (i + Then we wait for another migration 
from (i + to (i + 2,j) and so on so that in a finite time we reach (L + We 
work in the same way for the leftmost particle. We conclude that if <j>A = for all ei > 0, 
A>0, l<./V<oo there exists T ei such that £,t( x ) > Na for each x G [— 2L, 2L] 2 with 
probability larger than 1 — e%: hence 

P((0,0) is wet) > 1 - ex. (7.14) 

Suppose (pA > 0. For each e > there exists ei and T ei large so that (|7.14h holds and <pA 
small so that the probability of a death before T ei is as small as we want. Therefore 

P((0,0) is wet) > 1 - e. 

We conclude that for all L, e > and (m, n) G M the event G m ^ n = {(m, n) is wet}, which 
is measurable with respect to the graphical construction in -B m ,n, satisfies P(G m)Tl ) > 1 — e 
by taking T large and <j)A small. By comparison arguments with oriented percolation we 
get the result. 

ii) The idea is that even for 4>a small, there exists eft large so that the probability that the 
population size reaches N and then one individual migrates is small. One can prove the 
result by repeating the steps we did to prove Theorem 14.11 □ 
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7.3 Model III 

7.3.1 Proof of Proposition [5TTI 

We check the sufficient conditions for stochastic order from Theorem 12. II We call {P' ,V. .} 
the rates of (£t)t>o and {P.'.,F[ .} the ones of (r)t)t>o- They are given by 



PZ 



a<j)A,i if a < N A ~_ x _ \ a<f)A,2 if a < N A D i 



Pi = Pi = if < N 
a<p 2 if a > N A ; 1 p 



a<pi if a > Na; 
a,p ~ L a,p 

Let a < 7, ft < 5. We evaluate the terms in Condition (12.101) . The birth rates give 



R =Tt fl = A if a - A; > TV - M and /3 + £; < iV. 



E P P = %>5- ( 9+ii}- P y3 - /31{/3=5<Ar,ii=0} 
keX-.kyS-p+n 

E P 5 =1 {l>ji} P S = ^{j 1= 0,S<N} 



thus 



E n< E p l 

kex-.kys-p+j! iex-.l>j! 



(7.15) 



The death rates give 

E =1 {l>7-Q+^i}- f V 1 = 7%7=a,/ii=0}( < ^4,ll{7<AU} + 01 ^{N A <l} 

leX:l>"/-a+h L 



{a<AT A } + 02l{AT A <a} 



keX:k>hi 



thus 



E p a k > E ^ 



(7.16) 



k<EX:k>h± 

Now we consider the migration rates 

E = E ^^{fe<(a-Ar+M)A(W-y3)}; 

fee/a fceia 

By (|21))> - M> . setting l = k-8 + fi 



leX:l>"/-a+h 1 



E^. 5 ~~ E / ^ ]1 {K(7-A r +^)A(A r -^i- 



E / ^^{fc<(a-A''+-M')A(Ar- ( 3)} 
fceia 

A' 

=A |J{mi - 5 + £ > I > ji} n {0 < I < (a - N + M - d + /3) A (N - 5)} 



<A 



K 



|J{7 - a + mi > I > j { } n {0 < I < (7 - iV + M) A (iV - 5)} 



i=l 



: E ^ 1 {'<(7-A r +M)A(Af-<5)} 
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since 5 > (3 and 7 > a. Therefore 



In a similar way we note that 



(7.17) 



^ r a,/3 — ^ ]1 {fc<(Q-A r +A^)A(A r -/3)}; ^2^7,S ~ ^{l<(~f-N+M)A(N-S)} 

k£i d k£i d iei c iei c 



then, by setting k = I — 7 + a, the sum X^{i<(y-N+M)A(N-S)} is equal to 

K 

\ [Jim - J + a > k > hi} (1 {0 < k < {a - N + M) A (N - 5 - 7 + a)} 

K 

U{8-P + m i >k>h i }n{0<k<{a-N + M)A{N-P)} 



i=l 



< A 



Al{jfc<( a -Ar+Af)A(Ar-/3)} 



since iV — <5 — 7 + a<iV — /3. Hence 



(7.18) 



We get Condition (12TTUD by using (^T5]l and (TTTTD and Condition (I2TT1T) from (jTTlH]) and 
(fTTTHD . □ 

7.3.2 Proof of Theorem (5TTJ 

We follow the idea in \22\ proof of Theorem 2]. We assume d = 2. If d > 2 the proof 
works in a similar way. We take N, M such that N — M > Na- We fix x € Z d and we 
start from an initial configuration 770 (x) = N — M and 770(2) = for each z ^ x. We 
prove that starting from 770, after a finite time there is a migration of the largest flock of 
M (Na < M < N — Na) individuals into a site y ~ x which will give birth to iV — M 
individuals in the new site with large probability. 

For each x £ Z 2 we consider a modification (f?f )t>o constructed through the graphical 
representation in J x := [x — 1, x + l] 2 such that r)f(z) = for each z I x and t > 0: we 
take into account births, deaths and emigrations from x, births and deaths on each y ~ x, 
but we replace migrations of k individuals from y ~ x to x by the death of k individuals 
on y. For y ~ x, let 



£2 



{There exists T < 00 such that r§(y) = N-M\rfc(x) = N-M,rft(z) = 0,Vz ~ x}. 



Note that rjt > 77^ (it follows by construction from the graphical representation, since rf 
is built from 77; alternatively one can check Conditions (|2.10p -( j2.1ip by Remark I2.5D . In 
particular before T the process rjf(x) behaves as r)t(x) without immigration, and 77^(77) 
behaves as rjt(y). Therefore if E X)V occurs, rjr(y) > N — M. 

To make a comparison with an oriented percolation model, we follow |16j : between any 
two nearest neighbor sites x, y in Z 2 we draw a directed edge from x to y, denoted by [x, y): 
we say that one edge is open if E x y happens. This defines a locally dependent random 
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graph since E X y depends only on the graphical representation in I x . The probability of 
the directed edge [x, y) to be open is the same for all edges [x, y) and E XiV and E Z) t are 
independently open if x ^ z. 

We prove that for each e > there exists N large enough for W(E x y is open) > 1 — e. 
By translation invariance we suppose x = 0. We prove that the following events happen 
with large probability: first of all, starting from N — M, the number of visits to N — M + 1 
of r]f(x) := £t before visiting N A is at least iV 3 (Lemma [73]); if there are at least iV 3 
visits to N — M + 1, there are at least N 2 visits to N (Lemma I7.6P before reaching Na', 
if there are at least N 2 visits to N, there are at least iV 1 / 2 mass migrations of M > Na 
individuals to a fixed site y ~ (Lemma I7.7|) : finally one of these mass migrations gives 
birth to N — M individuals on y before reaching Na with large probability. 

I. First of all we prove that the number of visits R^ N M to N — M + 1 before reaching Na 
of the process (£t)t>o starting at N — M is large with large probability. 

Lemma 7.5 

lim F(R% M > iV 3 ) = 1. (7.19) 

N— >oo 



Proof . We construct a process ((t)t>0 with state space A := {Na, Na + 1, . . . N — M + 1} 
by coupling with in the following way: 

• jfN A <£t<N-M + l, then Ct = Cf, 

• i££t>N-M, then ( t = N-M + l; 

and Na is an absorbing state for (Ct)t>o- Each time that Ct hits N — M + 1, (an event 
which can happen only from below, that is if Ct moves from N — M to N — M + 1), so 
does Therefore we count the number of visits pS n m to N — M + 1 of the process Q 
starting at N — M + 1. Note that comes back to state N — M after visiting N — M + 1 
at an a.s. finite time Tm which satisfies 

H T M >t)< e~ xt (7.20) 

for each N, since if a mass migration of £t — (N — M) particles occurs then ^ comes back 
to N — M with rate A. The skeleton of the process (Ct)t>o moves as a discrete time random 
walk on A which comes back to N — M after visiting N — M + 1 with probability one, 
probability of birth p = 1/(1 + (ft) and probability of death 1 — p. We prove that 

lim P(i4 M > iV 3 ) = 1. (7.21) 

The probability that, starting at N — M, Q returns to N — M + 1 before visiting Na is 
given by Lemma 17.11 with r\ = Na, T2 = N — M + 1, j = N — M, q/p = (ft. Since after 
visiting N — M + 1 the walk returns to iV — M, by the Markov property (P/v_j\/+i (N — M) 
is the notation in Lemma l7.ip 

HR C N, M > N ") =(Pn-M+i(N - M)) N ' > (1 - <ftN-M-N A) N* 
>exp(-CN 3 (ft N - M ^) 

so that ([7^T]) (and then (fTTTU]) ) follows since (ft < 1. □ 
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II. Let R N be the number of visits of (£t)t>o to N before visiting Na starting at N — M. 
Lemma 7.6 

lim F{R% > N 2 ) = 1. (7.22) 

N— >oo 



Proof . By (l7~T9]l 



»(R N < N 2 ) =¥(R N < N 2 \R NM > N 3 ) +o(l) (7.23) 



where lim o(l) = 0. We define a family of i.i.d. random variables {Xj} i=1 ^3 such 

N-too '"' 

that Xi = 1 if & reaches N before TV - M at the i th visit to N - M + 1, otherwise. 
One possibility for X{ to be one is the birth of N individuals without any death or mass 
migrations. Such an event has probability larger than 

PN : ~ \i + 4> + xm/n) - \i + 4> + \m) ~ :p 

which does not depend on N. Therefore if Y is a binomial random variable with parameters 
p and N 3 , then P(X)fli^t < < W(Y < ^ 2 )> which converges to zero as N goes to 
infinity by the central limit theorem. □ 

III. Step II states that for each e > we are able to take iV large enough so that with 
probability larger than 1 — e the process £t reaches N at least N 2 times. We prove that in 
this case, for a fixed y ~ 0, with large probability there is a migration En = Ej^(0,y) of 
M individuals from to y at least N 1 / 2 times. 



Lemma 7.7 



lim F(E N > N 1/2 ) 

N^rOO 



Proof . Notice that when £t visits N there is a migration of M individuals from onto 
site y with rate X/(2d): if this is not the case, either a death at x or a different migration 
(that is less than M individuals onto y or a migration onto z ~ x, z 7^ y) occurs with rate 
smaller than N0 + \M(2d - l)/(2d) + (M — l)A/(2d). Thus the probability of a migration 
to y of M particles is larger than A/(2d(AM + N<f>)). 

The rest of the proof is identical to Step 2 of |22[ proof of Theorem 2]: the key point is 
that conditioning on {R^ N > N 2 }, En is larger than a binomial random variable V/v with 
parameters N 2 and X/(2d(XM + Nc/))), such that (V N - E(V A r))/(A^ 1 / 2+a ) converges to 
in probability for all a > 0. The claim follows by taking a € (0, 1/2). □ 

IV. We show that given at least N 1 / 2 emigrations from to y of M > Na particles, at 
least one of these flocks of individuals generates at least N — M + 1 individuals on y 
before reaching size Na- Every time there is a migration of M individuals to y, since 
M = M(Na) > Na, the process (77° (y))t>o 1S a birth and death chain with transitions 

Vtiv) -> Vt(v) + 1 at rate Vt(y)^{N A <^(y)<N-M+i} 
Vt(v) -> Vt(y) - 1 a t rate r]t(y)<P^{N A < v °(y)<N-M+i} 

Take the same chain on {Na, • • • , 00}. Since < 1, the chain is transient; therefore there 
is a positive probability q(<f>) that starting at M > Na the chain will go on to infinity. The 
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claim follows as in Step 3 of \22\ proof of Theorem 2] , since visits are enough for the 
probability to reach N — M + 1 at least one time to approach 1. 

We conclude that for each e > there exists N and T X}V large such that E xy occurs in a 
finite time T XtV with probability larger than 1 — e. 

V. Finally we conclude the comparison with the oriented percolation model on Z 2 . We 
say that percolation occurs if there exists an infinite path of directed open edges {(xo = 
0,xi) = ex,(xi,X2) = &2, ■ ■ ■ i ■ ■ •}; that is such that E XitXi+1 occurs for i = 0,1,.... 
Suppose ?7q(0) = N — M. If e\ is open, then r/^(0) reaches N, migrates to x\ and gives 
birth to N — M individuals on x\ before dying out. Then also e<i is open, therefore starting 
from iff 1 (x\) = N — M, it reaches N, migrates to x<i and gives birth to N — M individuals 
on 22 before dying out and so on: this is also true for the process r\ t > r]f for each x, 
therefore the existence of an infinite path in the percolation model implies the existence 
of an infinite path of individuals. 

We begin with one individual at x £ Z d . For each 4>a < oo, with positive probability ^(0) 
reaches N — M before in a finite time and we can start our construction. 
In order to prove that the existence of an infinite path in percolation model has positive 
probability if ¥(E XtV ) is large enough, one can follow |16t Theorem 3.2] and compare the 
process to a a site percolation model. Here we need d > 2, otherwise the construction 
does not work. The idea consists in making a comparison with an oriented site percolation 
model on the square lattice with both edges from a site open with a given probability 
7r, which can be taken as large as we want by taking N large. Since for such a model 
percolation occurs if tt is large enough, (|16j). there is survival with positive probability. 
If ?7o € Otv, then the upper invariant measure u, which exists by attractiveness, is not 
concentrated on the Dirac measure 5q and the claim follows. 

ii) The proof is similar to the one of Theorem 13.21 (Step (ii)), then we skip it. □ 

7.4 Model IV 

7.4.1 Proof of Lemma 16.11 

The process is a particular case of the reaction-diffusion process introduced in jTJ Section 
13.2]: by following the same notations, the reaction part of the formal generator (|2.ip is 

with q x (rj(x),7](x) + k) = rj(x)l^ = iy + 4>r](x)i^ = _iy The diffusion part is 

= : £ p(x,y)\c x ( V (x))lf(S^ y 1 r ] ) - f^)} 

x,y£Z d 

where p(x,y) = l{ y ^{2d)~ l and c x (r){x)) = %{ v (x)>N}- 

Since the maximal number of particles involved in a transition is finite and the birth and 
death rates grow linearly, the hypotheses of [71 Theorems 13.17 and 13.19] are satisfied 
hence existence and uniqueness of this process follow. □ 
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7.4.2 Proof of Theorem ETT1 



i) First of all we prove that there is stochastic order between Model / and Model IV. 
We consider the Model I as a process constructed on = Z^ d with birth rates null if the 
number of particles in a site is larger or equal to N. 

Lemma 7.8 Let ^ = $,t{4>' ■> A ') ^ e a process defined by {\2.1^ with rates given by {\°3. that 
is a Model I type process. Let rjt = n t (4>, <fi, A) be a Model IV type process. If <f> = $ , A = A' 
and £o(x) < N for each x € Z d ; then (j]t)t>0 is stochastically larger than (£t)t>o- 

Proof. Let (P;.,T..) and (P'.,r;.) be respectively the transition rates of (£t)t>o and 
{Vt)t>o- Note that an increase of particles in a site x with £t(x) = N is not possible, 
therefore £t( x ) — ^ for each x £ Z d and f > 0. 

We check conditions in Remark 12. 4[ Given < £(y) = r](y), 

P l(y) + r $(*U(iO = ^)%(</)<^-l} + A %(z)=7V,£(s/)<iV} ^ ^(y) + A1 {'?^)>A f ,'?(j/)<A f } 



-p 1 _i_ r 1 

77(y) ^ vmrty)' 



p m =*(f)%v)<"-i} < v(y) = Pr l{ 



and Conditions f|2.13|) f)2.14j) are satisfied. 

If £(x) = rj(x) (which is possible only if rj(x) < N), £(y) < r](y), 



P. 



\{x) + T l(x),av) + X1 {H(x)=N,t;(y)<N} > <Pv(x) + x1 {v(x)>nMv)<n} 

= p— i _i_ r 1 

Pg) =«(x)>0i 7 (x)=P^ )> 

so that Conditions (j27THjt — <[27THfl hold. □ 

Therefore by Theorem 13 . 21 there exists C (A, iV) such that if <p < (f> c (X, N) there is a positive 
probability of survival for Model I, and hence for Model IV. By taking eft > 1 one proves 
as in Model I (Step (ii) in proof of Theorem 13. 2p . that the process dies out: the existence 
of the critical parameter <p c follows from monotonicity with respect to <j>. 

ii) We skip it, since as in step i), stochastic order and Theorem 15.11 induce survival of the 
process. 

We prove that even if the process survives the expected value on each site is finite. Let 
77q (x) > N for each x € Z d and let (r)t)t>o be a process with N immortal particles per 
site, that is with transition rates 

T] t (x) —>rj t (x) + 1 at rate r] t (x) 

??f (a?) ->r] N (x) - 1 at rate ^nf (x)l {r] N {x)>N} . 

We define Ct(x) '■= ^(x) — N for each x G Z d , the birth and death process on N with 
birth rate N + (t(x) and death rate 4>(N + Ct(^))lfCt(aO>0}' Thus 

^E(C*(ar)) = E(Ct(x) + N) - 0E((Ct(x) + A0%(*)>o} < ^ " (£- l)E(Ct(*)) 
which implies 

E(C*(x)) <E(Co(x)) + iV/(0-l). 
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Therefore if Co(^) < n -> there exists c = c(n,N,cp) such that K(Q(x)) < c for each t > 
and x. The claim follows by taking C = c + N. □ 
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